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NOMENCLATURE 


[A] 

b 

C0 

Er , Eg , Ez 
FZ 

FZP 

FZQ 

FZT 


inner  radius  of  orthotropic  layer 
the  material  compliance  matrix 
outer  radius  of  orthotropic  layer 

ratio  of  inner  to  outer  radius  of  orthotropic  layer  (C0  =  a/b) 
engineering  moduli i  of  orthotropic  layer 

axial  force  on  orthotropic  layer  to  enforce  plane-strain  boundary 
coundi tions 

layer  axial  force  contribution  per  unit  of  internal  pressure 

layer  axial  force  contribution  per  unit  of  external  pressure 

total  axial  force  on  multiorthotropic-layered  cylinder  to  enforce 
plane-strain  boundary  conditions 

material-geometry  constant  from  hoop  strain  equivalence  condition 
(1  <  i  <  N- 1 ,  1  $  j  <  3) 


[J] 

[JI] 

k 

N 

P 

[Q] 

pl 

q 

r 

[R] 

RAP,  TAP, ZAP 

RAQ. TAQ, ZAQ 

RBP.TBP, ZBP 


matrix  of  G-jj  values  correctly  positioned 
inverse  of  [J] 

orthotropic  material  parameter 

total  number  of  layers  in  multilayered  cylinder 

internal  pressure  on  layer 

unknown  external  pressure  vector 

interference  pressure 

external  pressure  on  layer 

radius 

result  vector 

magnitude  of  r,9,z  stresses  evaluated  at  r=a,  caused  by  unit 
internal  pressure 

magnitude  of  r,8,z  stresses  evaluated  at  r=a,  caused  by  unit 
external  pressure 

magnitude  of  r,8,z  stresses  evaluated  at  r=b,  caused  by  unit 
internal  pressure 
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RBQ,TBQ,ZBQ  -  magnitude  of  r,0,z  stresses  evaluated  at  r=b,  caused  by  unit 
external  pressure 

Bij  -  compliance  for  cylindrical  problem 

6  -  radial  interference  of  two  multilayered  cylinders 

6-j  -  decrease  in  outer  radius  ' b '  of  inner  cylinder  per  unit  of 
external  pressure  (interference  loading  case) 

<50  -  increase  in  inner  radius  'a'  of  outer  cylinder  per  unit  of 
internal  pressure  (interference  loading  case) 

6U  -  radial  interference  reduction  per  unit  pressure  at  interface  of 
two  multiorthotropic-layered  cylinders 


e  -  strain 


vrd>v9z'vzr  -  Poisson's  ratio  of  orthotropic  layer  (e.g.  ,  vrg  = 
r-stress  producing  contraction  in  0-direction) 


e0 

^r 


a  -  stress 


[  ]  -  matrix  of  values 

(i)  -  pertaining  to  ith  orthotropic  layer 
Subscripts 

a,b  -  inner  and  outer  radial  evaluation  points 
r,9,z  -  radial,  hoop,  and  axial  directions 

1,2,3  -  directions  for  orthogonal  coordinate  system  (for  a  cylindrical 
system  1,2,3  correspond  to  r,0,z  directions) 


INTRODUCTION 


The  attractiveness  of  many  composite  materials  is  their  high  specific 
stiffness  making  them  an  ideal  choice  for  lightweight  applications.  In  order  to 
effectively  design  cylinders  that  incorporate  composites,  it  is  important  to  be 
able  to  predict  the  stress  distribution  caused  by  various  loadings  applied  to 
the  cylinder.  In  the  early  1960s,  Lehknitskii  (ref  1)  published  a  generalized 
plane-strain  stress  solution  for  a  monolayered  anisotropic  cylinder  unaer  com¬ 
bined  loadings  of  internal  pressure,  external  pressure,  and  axial  force.  More 
recently,  O'Hara  (ref  2)  simplified  these  equations  for  the  special  case  of  an 
orthotropic  material  and  cast  them  in  a  form  that  is  convenient  to  use.  In  real 
appi i cat  ions ,  however,  composite  cylinders  are  often  constructed  by  winding 
fibers  or  laying  up  fibers  at  various  angles.  Generally,  the  cylinder  is 
constructed  by  building  up  fiber  windings  in  positive  and  negative  wrap  angle 
pairs.  Each  of  these  positive  and  negative  wrap  angle  pairs  can  be  viewed  as  a 
single  orthotropic  layer.  The  whole  structure  can  be  considered  a  multi- 
orthotropic-layered  cylinder.  For  these  types  of  multilayered  cylinders,  it 
becomes  important  to  be  able  to  obtain  the  stress  aistribution  so  that  design 
and  analysis  can  be  pursued.  By  using  Lehknitski i 's  monolayered  solution  and 
the  proper  boundary  conditions,  the  multilayered  solution  can  be  constructed. 

The  equations  presented  herein  are  for  a  multiorthotropic-layered  cylinder  under 
plane-strain  boundary  conditions  with  internal  pressure,  external  pressure,  and 
interference  loadings. 

GEOMETRY 

A  multilayered  cylinder  can  be  viewed  as  an  assembly  of  many  single- layered 
cylinders.  It  is  fitting,  therefore,  to  begin  with  a  review  of  the  geometry  of 
the  monolayered  cylinder  problem. 
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For  the  raonolayered  case,  the  cylinder  is  assumed  to  be  long  with  ends  that 
are  fixed  (see  Figure  1).  The  fixed-end  condition  implies  zero  axial  strain 
through  the  radial  thickness  of  the  cylinder.  The  cylinder  has  an  inner  radius 
'a'  and  an  outer  radius  'b'.  A  cylindrically  orthotropic  material  is  assumed 
with  its  principal  axes  coincident  with  the  cylinder.  Orthotropic  material,  in 
general,  is  characterized  by  nine  independent  material  constants.  These 
constants  consist  of  three  engineering  modulii  (E},E2»E3)/  three  shear  modulii 
(g12 '^23 < G31 ) '  and  three  Poisson's  ratios  (Pi2»y23*w31 ) •  The  numbers  1,2,3 
indicate  the  principal  material  directions.  For  the  above  assumptions,  the 
1,2,3  directions  correspond  to  the  radial,  hoop,  and  axial  directions  of  the 
cylinder  (r,9,z).  The  cylinder  can  be  subjected  to  internal  pressure  'p'  and 
external  pressure  ’ q ’ .  In  addition,  since  the  principal  stress  directions 
correspond  to  the  principal  geometry  directions  of  both  the  cylinder  and  the 
applied  loadings,  shear  effects  are  eliminated. 


STRESS  EQUATIONS  FOR  MONOLAYERED  ORTHOTROPIC  CYLINDER 


Lekhni tski i ' s  solution  (ref  1)  for  a  cylinder  with  one  anisotropic  layer  can 
be  simplified  for  the  case  when  the  material  is  orthotropic.  This  simplifica¬ 
tion  was  done  by  O'Hara  (ref  2)  and  resulted  in  three  stress  equations  that 
correspond  to  the  r,9,z  directions  and  one  equation  that  predicts  the  axial 
force  necessary  to  enforce  the  fixed-end  constraint.  These  equations  are  given 
below  and  are  identical  to  those  found  in  Reference  2  with  one  correction  for  a 
typographical  error  in  the  az  equation.  Also,  the  axial  force  name  has  been 
changed  from  PP  to  FZ  for  clarity. 
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where  ’ r '  is  the  radial  position  in  the  cylinder, 


(4) 


C0  = 


a 

b 


(5) 


and  the  components  A^ j  are  the  elements  of  the  compliance  matrix  as  given  in 
Hooke's  Law, 


[e]  =  [A]  [a] 

' k '  is  an  orthotropic  material  constant  given  by 


where,  in  general. 


^ij  ■  Ai j 


5ll 

322 


'33 
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In  developing  the  multilayered  solution,  applying  the  correct  boundary 
conditions  at  the  interface  of  two  orthotropic  layers  necessitates  the  use  of 
stress  values  evaluated  at  the  inner  (r=a)  and  outer  (r=b)  surfaces  of  each 
layer.  This  evaluation  process  leads  to  the  six  stress  equations  given  below. 
Stress  Equations  Evaluated  at  Inner  Radius  (r*a) 

ar<a  =  RAP-p  +  RAQ-q  (9) 

RAP  =  -1  ,  RAQ  =  0 

ae,a  =  TAP-P  +  TAQ-q  (10) 


Ml+C^kj 

(1-C02k) 


TAQ  = 


-2kC0k-l 

(1-C02k) 


az,a  =  ZAP-p  +  ZAQ-q 


ZAP  * 


U13-A23  -TAP) 
*33 


ZAQ  = 


:A23;TAQ 

a33 


Stress  Equations  Evaluated  at  Outer  Radius  ( r=b) 

ar,t)  ~  RBP-P  +  RBQ-q 
RBP  =  0  ,  RBQ  =  -1 

<?0  ,b  =  TBP'P  +  TBQ '  <=1 


TBP  =  -TAQ • C02  ,  TBQ  =  -TAP 


(ID 


(12) 


(13) 


aZii)  a  ZBP'P  +  ZBQ-q 


ZBP  =  -ZAQ*C08 


ZBQ  =  -ZAP  + 
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^33 


Also,  the  axial  force  equation  can  be  rewritten  in  a  similiar  form. 


(14) 
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FZ  =  FZP-p  +  FZQ-q 


(15) 


2rra£ 

[ ( TAP+TAQ+1 ) (A23_A^3 ) ] 

a33 

(1-k*) 

a23 


FZQ  = 


2nb  l 

[ (TBP+TBQ+1 ) (A23-A13 ) ] 

a33 
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'23 


In  each  of  the  six  stress  equations  above,  there  are  two  constants  (e.g., 
TAP,  TAQ )  which  determine  the  magnitude  of  the  contribution  to  the  given  stress 
caused  by  the  internal  and  external  pressure.  These  constants  are  material- 
and  geometry-dependent,  and  each  has  a  three-letter  name.  The  first  letter 
corresponds  to  the  stress  direction  (r,0,z),  the  second  letter  signifies  the 
ooint  at  which  the  stress  has  been  evaluated  (a,b),  and  the  third  letter 
corresponds  to  the  pressure  it  applies  (q,p).  For  example,  TAQ  is  the  constant 
for  the  Theta  stress,  evaluated  at  r  =  a  (A),  magnifying  the  external  pressure  q 


(9). 


STRESS  SOLUTION  FOR  MULTIORTHOTROPIC-LAYERED  CYLINDER 

In  this  section  the  methodology  of  constructing  the  stress  solution  for  a 
multilayered  cylinder  (see  Figure  2)  is  discussed  where  each  layer  has  material 
and  geometry  definitions  identical  to  those  discussed  for  the  monolayered  case. 

The  first  steo  to  construct  the  multilayered  solution  is  to  equate  the 
circumferential  (hoop)  strain  of  adjacent  orthotropic  layers  at  their  interface. 
In  other  words,  we  equate  the  hoop  strain  at  r=b  of  the  i^*1  layer  with  the  hoop 
strain  at  r=a  of  the  i+1  layer 

e9,b(i)  =  e9,a(  i+D  (16a) 
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Recalling  that  [e]  =  [A]  [a]  and  that  1,2,3  correspond  to  the  r,9,z  directions 

e2,b(i)  =  e2,a(1+D  (16b) 

or 

A12(i )  ■<7r,b('i )  +  a22(‘>)  •<7e,b('i)  +  A23*az,b(  "> )  = 

A12(i  +  1)  ,<7r,a('i  +  l)  +  A22('i+1),a0,a('i+1)  +  A23‘az,a( i  +  1 )  (I7) 

Now,  substituting  the  expressions  for  stress  found  in  Eqs.  (9)  through  (14) 
we  have 

-A12(i)-d(i)  +  A22(i)-[TBQ(i)-q(i)  +  TBP(i)-p(i)]  + 
A23(i).[ZBQ(i)-q(i)  +  ZBP ( i )  •  p ( i )  ]  = 

-A12 ( i +1 ) • p ( i+1 )  +  A22(i  +  1) • [TAQ( i  +  1) -q( i  +  1 )  +  TAP ( i  +  1 ) *  p( i  +  1 )  ]  + 

A23 ( i+1 ) • [ZAQ( i+1 ) *q( i+1 )  +  ZAP( i+1 ) *p( i+1 ) ]  (18) 

Rearranging  and  noting  from  equilibrium  conditions  that 

P(i)  =  q(i-l)  (19) 

we  have 

Gil*q(i-1)  +  G-j2'c)('i)  +  G  ^  3  •  q  ( i  +- 1 )  =  0  (20) 

where 

Gil  =  -C0=  ( i }  */322(  ■» )  *TaQ(  i ) 

G-j2  =  -[/3i2('')  “  /3i2(i  +  l)  +  ^22(  ^ )  "TAP(  i )  +  022  ( i+l )  ' TA^  ( +  D  ] 

Gi3  =  -j322  ( i+l ) 'TAQ( 'i  +  l ) 

For  each  two-layer  combination  there  is  one  equation  defining  the  hoop  strain 
equivalence  condition.  For  a  cylinder  with  ' N ’  orthotropic  layers  there  are 
N-l  equations  and  N-l  unknowns.  Setting  up  these  equations  in  matrix  form 
gives  the  following: 
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G11  g12  g13 

“q(0)‘ 

g21  g22  g23  0 

q  ( l ) 

G31  g32  g33 

q(2) 

0  - 

. 

gN-1,1  gN-1,2  gN-1,3 

q(N) 

0 


(21) 


Since  q ( 0 )  and  q(N)  are  simply  the  internal  and  external  pressure  of  the  overall 
cylinder  and  are  known  quantities,  the  above  matrix  expression  reduces  to 


G12  g13 

q  ( 1 ) 

’  -Gn-p(i) 

g21  g22  g23  0 

q(  2 ) 

0 

g31  g32  g33 

q  (3 ) 

, 

s 

0 

0  GN-2 , 1  gN-2 , 2  gN-2 , 3 

GN-1,1  gN-1 , 2 

q(N-l) 

—  — 

_GN-1 , 3 -q(N) 

(22) 


[j]  CQ]  =  [R] 

Solving  for  the  matrix  [Q]  we  have 


(23) 


[Q3  =  [J]”1  [R]  =  [JIJ  [R] 

where 


(24) 


[JIJ  =  [JJ-1 


(25) 
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Once  [JI]  is  determined,  the  interface  pressures  can  be  found  easily  by  the 
following  equation: 

For  (1  $  i  <  N-l): 

q(i)  =  -Jl-j ,l'Gii*p(l)  -  JIi ,N-l*GN-l,3'q(N)  (26) 

By  recalling  that  p(i)  =  qp-1),  we  can  now  calculate  the  stress  distribu¬ 
tion  in  each  layer  by  using  the  general  stress  equations  (Eqs.  (1)  through  (4)) 
with  p ( i )  and  q ( i )  as  input. 

In  addition  to  finding  the  stress  distribution  in  each  layer,  the  axial 
force  necessary  to  constrain  each  layer  to  zero  axial  strain  can  be  calculated 
using  Eq.  (15)  with  p ( i )  and  q ( i )  as  input.  By  summing  the  axial  forces  on  each 
layer,  we  can  determine  the  total  axial  force,  FZT,  on  the  multilayered  cylinder 

N 

FZT  =  ^  FZ( i )  (27) 

i=l 

MULTILAYERED  SOLUTION  FOR  PRESS  AND  SHRINK  FIT 

For  many  engineering  applications,  it  is  desirable  to  assemble  two  cylin¬ 
ders  by  shrinking  or  press  fitting  one  part  upon  the  other.  This  process  pro¬ 
duces  a  pressure  at  the  mating  surface  of  the  two  vessels.  The  stresses 
resulting  from  this  pressure  on  each  cylinder  can  easily  be  found  using  the 
method  discussed  in  the  previous  section.  To  produce  this  mating  pressure,  the 
male  member  must  have  an  outer  radius  greater  than  the  inner  radius  of  the 
female  member.  The  difference  in  these  two  radii  is  called  the  interference  6, 
and  is  the  radial  deformation  experienced  by  both  cylinders  at  the  interface. 

If  a  unit  pressure  is  applied  to  the  outer  diameter  of  the  male  or  inner 
cylinder,  a  decrease  in  the  outer  radius  of  6 -j  results.  Likewise,  a  unit 
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pressure  applied  to  the  inside  of  the  female  or  outer  member  produces  an 
increase  in  its  inner  radius  of  50.  Therefore,  the  reduction  in  interference 
per  unit  of  interface  pressure  is  given  by 

5U  =  «i  ♦  «o  (28) 

If  the  interference  is  at  the  i,  i+1  layer  interface, 

6U  =  -b(  i )  •£0,b('i )  +  aM  +  1)  *€0,a(i  +  l)  (29) 

In  general,  however,  the  interference  is  small  compared  to  b(i)  and  a(i+l), 
therefore 

b( i )  =  a( i+1 )  (30) 

and 

6U  =  b(i)-[eefa(i+l)  -  €0,b(i)]  (31) 

The  strains  £95(1)  and  £0(a(''  +  l)  are  easily  found  from  the  procedure 
discussed  in  the  previous  section.  Once  this  is  done,  the  interference  per  unit 
of  internal  pressure  must  be  scaled  by  pressure  Pj  to  equal  the  total  inter¬ 
ference 

6  =  Pl-5U  (32 ) 

Thus,  Pj  is  given  by 


Pj  is  the  interface  pressure  that  must  exist  to  produce  a  radial  defor¬ 
mation  at  the  interface  in  both  cylinders  equal  to  the  interference  6. 

RESULTS 

The  multilayered  solution  was  used  to  investigate  the  stress  and  strain 
distribution  within  a  ten-layered  cylinder  (N=10).  The  cylinder  had  an  inside 
radius  a(l)  ■  1  inch  and  an  outside  radius  b(10)  =  2  inches,  and  each  of  the  ten 
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layers  was  0.1  inch  thick.  The  material  used  was  an  IM6/epoxy  with  a  55  per¬ 
cent  fiber-volume  ratio.  The  layup  considered  consisted  of  five  hoop-axial 
pairs  starting  at  the  inside  radius  of  the  cylinder.  The  material  properties 
for  the  hoop  and  axial  fiber  orientations  are  given  in  Table  I.  The  results 
generated  consist  of  three  sets  of  plots  (Figures  3,  4,  and  5).  Each  set  of 
plots  corresponds  to  one  of  three  loading  conditions.  For  each  of  these  loading 
conditions,  four  graphs  were  generated  displaying  the  stress  and  strain  distri¬ 
bution  within  the  cylinder  as  a  function  of  radial  position.  These  graphs 
include  stresses  in  the  r,9,z  directions  and  strains  in  the  r,9  directions.  The 
three  loading  conditions  were  unit  internal  pressure  { p ( 1 )  =  1.0),  unit  external 
pressure  (q(10)  =  1.0),  and  a  1-microinch  interference  at  the  5  to  6  layer 
interface. 


TABLE 

I .  MATERIAL 

PROPERTIES 

FOR  IM6/EP0XY 

55%  FIBER-VOLUME 

RATIO 

Fiber 

Direction 

Er  (Mpsi) 

Eg  (Mpsi) 

Ez  (Mpsi) 

Vre 

u9z 

vzr 

Hoop 

1.126 

23.31 

1.126 

0.0152 

0.3147 

0.3991 

Axial 

1.126 

1.126 

23.31 

0.3991 

0.0152 

0.3147 

The  multilayered  solution  was  also  compared  to  a  finite  element  solution 
for  the  cylinder  with  the  first  two  loading  conditions  mentioned  above.  The 
ABAQUS  finite  element  code  was  used  to  produce  stress  results  used  in  the  com¬ 
parison.  The  finite  element  model  contained  20  eight-node  quadratic  axisym- 
metric  elements  with  two  equal-sized  elements  used  to  model  each  of  the  ten 
orthotropic  layers.  The  comparison  was  limited  to  the  radial,  hoop,  and  axial 
stress  values  evaluated  at  the  inner  and  outer  radii  of  each  orthotropic  layer. 
Tables  II  and  III  contain  the  results  of  this  comparison. 
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COMPARISON  BETWEEN  MULTILAYERED  SOLUTION  AND  ABAQUS  FINITE  ELEMENT  SOLUTION  USING  20  CAX8  ELEMENTS 
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TABLE  III.  COMPARISON  BETWEEN  MULTILAYERED  SOLUTION  AND  ABAQUS  FINITE  ELEMENT  SOLUTION  USING  20  CAX8  ELEMENTS 
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DISCUSSION 


The  stress  and  strain  distributions  for  the  three  loading  conditions  are 
displayed  in  Figures  3,  4,  and  5.  In  each  case,  the  radial  stress  must  be  equal 
to  the  negative  of  the  applied  pressure  at  the  loading  surface.  Therefore,  for 
case  1  involving  internal  pressure,  the  radial  stress  was  equal  to  -1.0  psi  at 
the  inside  radius  and  zero  at  the  outside  radius.  For  case  2  involving  external 
pressure,  the  radial  stress  was  equal  to  -1.0  psi  at  the  outside  radius  and  zero 
at  the  inside  radius.  And  for  case  3,  the  radial  stress  was  equal  to  zero  at 
the  inner  and  outer  surfaces  and  equal  to  about  -0.98  DSi  at  the  5  to  6  layer 
interface,  -0.98  was  the  interference  pressure  associated  with  the  1-microinch 
i nterf erence .  In  the  three  graphs  showing  the  hoop  stress  distribution  for  the 
three  loading  conditions,  it  is  clearly  seen  that  the  hoop  layers  took  up  most 
of  the  load.  This  is  primarily  due  to  the  low  radial  stiffness  of  botn  the  hoop 
and  axial  layers  (ref  2).  It  is  also  observed  that  the  axial  stress  for  the 
three  loading  conditions  was  relatively  small  when  compared  to  the  hoop  stress. 
This  was  expected  because  the  axial  stress  was  a  second  order  effect  related  to 
the  Poisson's  contraction  and  axial  stiffness  of  each  layer.  As  a  result  of  the 
combination  of  the  low  axial  and  radial  stiffness  of  the  hoop  layer  and  the  low 
radial  stiffness  of  the  axial  layer,  a  very  small  axial  load  was  needed  to 
enforce  the  plane-strain  boundary  condition. 

The  graphs  of  the  hoop  strains  for  the  three  loading  conditions  show  some 
interesting  results.  The  hoop  strain  for  the  internal  pressure  loading  was 
smooth  and  monotonically  decreasing,  while  the  hoop  strain  for  the  external 
pressure  loading  was  smooth,  but  somewhat  larger  in  magnitude  a  small  distance 
away  from  the  inner  surface  than  at  the  inner  surface.  The  hoop  strain  for  the 
interference  loading  case  shows  a  discontinuity  at  the  interface  of  loading. 
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This  was  expected  because  the  interference  pressure  causes  a  shrinking  of  the 
outer  diameter  of  the  inner  cylinder  and  an  expansion  of  the  inner  diameter  of 
the  outer  cylinder.  Finally,  the  radial  strains  for  the  three  loading  cases 
were  not,  in  general,  continuous.  The  radial  strain  for  the  interference 
loading  shows  this  most  clearly  at  a  radial  position  of  1.6  inches.  This  is 
because  the  radial  strain  depended  on  both  the  material  properties  and  stresses 
of  each  layer  at  the  interface  of  each  layer. 

The  theoretical  and  finite  element  comparisons  of  the  radial,  hoop,  and 
axial  stresses  at  the  inner  and  outer  surfaces  of  each  orthotropic  layer  are 
shown  m  Tables  II  and  III.  Table  II  shows  the  comparison  for  the  unit  internal 
pressure  loading  case,  and  Table  III  shows  the  comparison  for  the  unit  external 
pressure  loading  case.  For  both  loading  cases,  the  comparison  shows  excellent 
agreement  between  the  theory  and  the  finite  element  results.  For  the  finite 
element  cases,  the  stresses  were  obtained  at  the  nodes  and  are  somewhat  less 
accurate  than  stresses  obtained  at  the  integration  points  of  the  element. 

CONCLUSION 

A  plane-strain  elastic  stress  solution  for  a  multiorthotropic-layered 
cylinder  has  been  developed.  The  solution  has  been  shown  to  be  in  excellent 
agreement  with  finite  element  results.  The  solution  has  also  been  implemented 
in  a  FORTRAN  program  to  be  used  as  a  tool  for  the  designer. 
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Theoretical  stress  and  strain  plots  for  N  =  10,  5*(l-hoop,  1-axial) 
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Theoretical  stress  and  strain  plots  for  N  =  10,  5*(l-hoop,  1-axial) 
(see  Table  I),  a(l)  =  1.0  in.,  =  2.0  in.,  1.0  pin.  interference 
at  the  5  to  6  layer  interface. 
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